Characteristic length scale and energy of a vortex line in a dilute, superfiuid Fermi gas 
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We calculate the characteristic length scale of a vortex line in a dilute, superfiuid gas of fermionic 
atoms, and find that it is in general smaller, and has a weaker density dependence than the BCS 
coherence length. Taking this into account, we find the energy of a vortex line to be larger than 
has previously been estimated. As a consequence, the critical frequency for formation of vortices in 
trapped fermion gases can be significantly larger than earlier calculations have suggested. 



The transition to a superfiuid Bardeen-Cooper- 
Schrieffer (BCS) state at low temperatures is a generic 
feature of Fermi systems with an attractive effective in- 
teraction. The low-temperature superconductors and liq- 
uid 3 He are the most familiar examples, but the phe- 
nomenon can also be seen in more exotic systems like 
atomic nuclei, and superfluids are also thought to be im- 
portant in neutron stars [0. Shortly after Bose- Einstein 
condensation in a dilute gas of 87 Rb atoms was achieved 
the possibility of observing the BCS transition to a 
superfiuid state in a dilute gas of trapped fermionic atoms 
was proposed J3|. Trapping and cooling of fermionic al- 
kali atoms have now been achieved, reaching tempera- 
tures as low as - T F /4 for 40 K @,| and 6 Li §-§, where 
Tp is the Fermi temperature. The typical transition tem- 
perature for weakly interacting fermions is T c ~ 10~ 3 Tf, 
but one hopes to create systems with T c ~ 10 _1 Tp in the 
experiments. If temperatures below T c can be reached, 
a new and exciting laboratory for studying properties 
of fermion superfluids will be at our disposal. Finding 
experimental signatures of the BCS transition in these 
systems is therefore important, and considerable effort 
has been invested in this problem ||. 

The study of quantized vortices in Bosc-Einstein con- 
densates has led to many interesting results (To) , and re- 
cently the vortex state in trapped fermionic gases has also 
been studied Jll|,|l2| as a means for detecting a superfiuid 
state in these systems. A fundamental quantity in these 
considerations is the critical rotational frequency for the 
creation of a vortex line with one quantum of circulation. 
The gases in the experiments are inhomogeneous because 
of the trapping potential, but if the vortex is produced 
near the center of the trap, one can estimate the critical 
frequency starting from the energy per unit length, £ v , of 
a vortex in homogeneous matter Jl^ , p^[ . Recently, Bruun 
and Viverit [fl2| calculated £ v at zero temperature, as- 
suming that the size of the vortex core, the region where 
the order parameter drops from its value in homogeneous 
matter to zero, is given by the BCS coherence length £q- 
However, several studies of vortices in low-temperature 
superconductors |L4|-[l6| have found that the formation 
of bound quasiparticle states localized near the center of 



the vortex can cause the size of the vortex core to be 
much smaller than £o at low temperatures. The same 
effect has been found for vortices in superfiuid neutron 
matter , and we will in this Letter argue that it is 

also present in dilute gases, and that this can lead to a 
significantly larger vortex energy than that found in [fl2f . 

We consider a dilute, homogeneous gas of fermionic 
atoms in two different hyperfine states, held at temper- 
ature T = 0. In this limit the Pauli principle dictates 
that the effective s-wave interaction between two atoms 
in the same hyperfine state vanish. The interaction be- 
tween two atoms in different states can be approximated 
by the s-wave scattering length a. For a negative scatter- 
ing length, the interaction is attractive and if the num- 
ber of particles in the two states is the same, the T = 
ground state of the gas is superfiuid. Theoretically, the 
formation of a superfiuid state is signaled by a non-zero 
value for the energy gap Ao, which is proportional to the 
critical temperature T c . In the dilute gas limit k F \a\ <C 1, 
where k F is the Fermi wavenumber, it is given by 
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where e F — fi 2 k F /2m a , m a being the mass of an atom, 
is the Fermi energy [[l9 20 1. Corrections from the so- 
called induced interaction modify the numerical prefactor 
in this result, but the dependence on k F and \a\ remains 
the same f0|[20|. 

The energy per unit length of a vortex line with re- 
spect to the homogeneous ground state of the superfiuid 
can be estimated by adding the kinetic energy associated 
with the velocity field around the vortex, and the loss 
of condensation energy in the vortex core p"2| . In this 
approximation the vortex core is treated as a cylindrical 
column of normal matter with radius equal to the BCS 
coherence length £o = ^wf/ttAo, where v F = hk F /m a - 
The result is 
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where n a = hp/ 6tt 2 is the number density of one hy- 
perfine state, and D — 1.36. The energy can also be 
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calculated from Ginzburg-Landau theory, and the result 
turns out to be of the same form, but with D = 1.65 [fL2| , 
again assuming that £o is the characteristic length scale 
for the vortex. However, at low temperatures the micro- 
scopic properties of the system are important, and, as 
we will argue below, the existence of quasiparticle states 
localized near the center of the vortex has to be taken 
into account at very low temperatures. 

The microscopic properties of vortices in fermion su- 
perfluids can be derived from the Bogoliubov-de Genncs 
(BdG) equations , which we write as 

W Bd G*i(r) - £i*i(r), (3) 

where 

_ / H (r) A(r) \ 

n BdG - (^ A * (r) _ h:(y) J, (4) 

with H a = —h 2 \7 2 /2m a — ep, and = (iii,Vi) T , where 
U{ and Vi are the coefficients in the Bogoliubov quasi- 
particle transformation, and i symbolizes the quantum 
numbers characterizing the quasiparticle states. The or- 
der parameter A(r) at temperature T is determined by 
the self-consistency condition 

i(r)= i^ E „, (rK(r)tani (_!_), (5) 

where ks is Boltzmann's constant. In cylindrical coor- 
dinates (p,(p,z), a vortex with one quantum of circu- 
lation can be described by a complex order parameter 
A(r) = |A(r)|e ,e ( r ). For a straight vortex line along 
the axis of the container, |A(r)| = A(p), and 9 = —Lp, 
where A(p) starts out at zero in the center of the vortex 
core at p = 0, and increases asymptotically to the value 
in homogeneous matter Ao- The low- lying quasiparti- 
cle states can be obtained analytically, as first shown by 
Caroli et al. |2l| . Kramer and Pesch [Q used a simi- 
lar approach in a study of vortices in type II supercon- 
ductors, where they showed that at low temperatures, a 
new length scale, different from characterizes the size 
of the vortex core. Their calculation, translated to our 
notation, is as follows: We label the lowest eigenstates 
of the BdG Hamiltonian by their angular momentum p, 
equal to half an odd integer, and the z component of the 
momentum k z = kp sin 9, where 9 is the angle between 
kj? and the xy plane. For < p <C kptio, where kp£o 
is of order 10-100 in the weak coupling regime, and for 
small p, one finds that these eigenstates are of the form 

BIB 

?v(r) = (S^) ^ ip+1/2W J, + iMk F pcos9) (6) 
v,(r) = e-^-^J^ 1/2 (k F pco S 9) (7) 



where J p is a cylindrical Bessel function and L z is the 
length of the vortex line. We define the characteristic 
length scale of the vortex core via the behavior of the 
order parameter near p = 0: 

lim A(r) = e^Ao^-. (8) 

As T — > 0, the sum in Eq. (||) is dominated by the lowest 
positive energy state, p = 1/2. Inserting Eqs. (||) and 
(@) for p = 1/2, and using J p (x) ps ^ (f ) P for x < 1, we 
obtain the behavior of the order parameter near p = 
from Eq. (||) as 

A(r) = ^-^A e->^cos^anh(^) . 
m a L z v F X \ 2k B T ) 

(9) 

By comparing with Eq. (|^) we find 

- = cos^tanh ( , (10) 

where we have used Dp — kp/m a . Converting the sum 
over k z to an integral over 9 through the relation 

£- jT*«». (") 

we find 

^mf^w^jm, (12) 

since tanh(i? 1 / 2 .fc i /2fc_BT) ps 1 as T — ► 0. The character- 
istic length scale at low temperatures T <C T c is therefore 
given by 



and is seen to be a result of the dominating role played by 
the low-lying eigenstates of the BdG Hamiltonian at low 
temperatures. Numerical results obtained by Gygi and 
Schliiter JT^,[ll| have supported this conclusion. Further- 
more, De Blasio and Elgar0y |]l7| , ^8| , in a study of vor- 
tices in superfiuid neutron matter at T = 0, considered 
various ways of defining the vortex core size. They all 
gave characteristic sizes generally smaller than £q, the 
differences being particularly pronounced for low values 
of Ao, which is the case in the dilute gas limit. The 
most important feature, however, was that they had a 
much weaker dependence on kp\a\ than the BCS coher- 
ence length, and that the overall trend of the numerical 
results was well approximated by the Kramer-Pesch esti- 
mate for the vortex core size, £i. The ratio of this quan- 
tity to f , Ci/Co ~ e-*/ 2kF W/k F \a\, is plotted in Fig. 
and can be seen to be much smaller than one in the weak 
coupling regime. 
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however, in contrast with Eq. (g), D is now a function of 
through the ratio £i/£o- Since £1 <C £o hi the weak 
coupling regime, we have 
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FIG. 1. The ratio £i/£o of the Kramer-Pesch estimate for 
the vortex core size to the BCS coherence length. 



To estimate the vortex energy per unit length at zero 
temperature, we follow the approach of fj^ ], but with the 
important modification that we take the vortex core size 
to be given by £i instead of £o- The vortex is modeled as 
a cylindrical region of radius £i containing normal mat- 
ter, surrounded by a velocity field with magnitude given 
by v s = h/2m a p. The superfluid is confined within a 
cylinder of radius R c 3> £i, and is assumed to be uni- 
form along the z-direction, which also defines the vortex 
axis. The energy per unit length £ v of the vortex state 
with respect to the homogeneous superfluid can then be 
divided into two contributions. The first is the kinetic 
energy associated with the flow of the superfluid outside 
the vortex core, given by 
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and the second contribution, coming from the loss of con- 
densation energy in the vortex core, can be estimated by 
multiplying the BCS result for the condensation energy 
per volume, e con d = 3AQn CT /4e_F, with the area of the 
vortex core tt^ 2 . This gives 



£"cond — 



ith 2 n a 3 /£i 
2m a 7T 2 \£ 



(15) 



The total energy can therefore be written in the usual 
form 
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and the vortex energy is to a good approximation given 
by ln(i? c /£i), which has also been found to be the case 
in neutron matter jll] . Note that D is essentially the in- 
verse of the quantity plotted in Fig. ^, and so it increases 
rapidly with decreasing /c^|a|. At the limit of the region 
where weak coupling can reasonably be expected to hold, 
kp\a\ ~ 0.4, D « 20.0, which is an order of magnitude 
larger than in the estimate in Eq. (g) . 

We can now estimate the critical rotation frequency for 
a dilute gas of atoms trapped in a cylindrically symmetric 
harmonic oscillator potential 

V(r)^^m a u; 2 z [z 2 + X 2 T (x 2 +y 2 )}, (18) 

where At describes the anisotropy of the trap. We will 
use the Thomas-Fermi result for the density profile of 
this gas, given by 



n„(p,z) = n.o-,0 1 - 



X 2 T p 2 + z 2 ^ 3/2 

m 



(19) 



where ?v o = n CT (0, 0) is the central density of the cloud, 
R z = (48A f (T A|) 1 / 6 / osc is the extent of the cloud in the 
z-direction, N a is the number of atoms in the hyperfine 
state <7, and l osc = ^Jh/m a uj z J2j|. Following Lundh et 
al. jl3| , we divide the cloud into vertical slices of height 
dz, and use (|l6| ) for a cylinder of radius p\ such that 
£i <C pi <C R± = Rz/^t, where the gas can be consid- 
ered uniform. The energy per unit length of the slice at 
z can then be written as 
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with Rx(z) = (1 - z 2 /R 2 z ) 1 / 2 R z /\ T being the extent of 
the cloud in the p-direction for a given z, and n a (0, z) 
is the density along the z-axis. Inserting Eq. ( |l9] ) and 
integrating, one finds 
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where £i(0) is the value of £i at the center of the cloud. 
Integrating over z, we obtain the total energy of the vor- 
tex as 
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The critical rotation frequency lj c \ for creating a vortex 
in the trap is given by io c \ — E v /L v where L v = N a h 
is the angular momentum of the vortex state. Using 
N a = 7T 2 i?^n CTi0 /8Ay, we find 



3 I 2 ( R 
uj c1 = -u±-±-ln 0.379—!- 

2 ± R 2 ± \ a(o) 



(23) 



where l± = l OS c/^T- To compare with the result of Ref. 
JL2| , we take At = 1 so that uj z = uj±_ = oj, and set 
fc_F,o|a| = 0.4, ep = 200^, corresponding to an isotropic 
trap with N a ~ 1.3 • 10 6 . With these values, we obtain 
uj c \ w 0.014u, four times larger than in |l2| . In Fig. (|^) 
the ratio between the critical frequency found in this pa- 
per and the one found in jl2| is plotted as a function of 
kp\a\, with the other parameters kept at the same values 
as above. Note that lo c \ turns negative for small val- 
ues of kp\a\, reflecting that the condition £i <C R± is 
violated. For the parameters used in this example, this 
occurs at kp\a\ ~ 0.01 for u c i given by Eq. (|23|), and for 
kp\a\ « 0.33 for the corresponding result in Ref. u2]. 
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FIG. 2. The critical frequency, Eq. (|23J) (solid line), the 
critical frequency found in Ref. [10] (dashed line), and their 
ratio (dot-dashed line), for an isotropic trap with the param- 
eters given in the text. 



To conclude, we have estimated the energy of a vortex 
in a uniform, dilute gas of atoms at T = 0. The quasi- 
particle states in the core of the vortex introduce a new 
characteristic length scale £i which has a weaker den- 
sity dependence than the BCS coherence length £o and is 
much smaller than £o in the dilute gas regime kp\a\ 1. 



Taking this fact into account leads us to predict a con- 
siderably higher vortex energy than what was found in 
earlier calculations. For the case of a trapped gas of 
fermionic atoms, this means that the critical frequency 
for the formation of a vortex state can be 5 — 10 times 
higher than previously estimated. 
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